Multiband superconductivity in NbSe2 from heat transport 
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Abstract 

The thermal conductivity of the layered s-wave superconductor NbSe2 was measured down to T c /100 throughout 
the vortex state. With increasing field, we identify two regimes: one with localized states at fields very near H c \ 
and one with highly delocalized quasiparticle excitations at higher fields. The two associated length scales are 
most naturally explained as multi-band superconductivity, with distinct small and large superconducting gaps on 
different sheets of the Fermi surface. 
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In the classical theory of superconductivity all elec- 
trons on the Fermi surface contribute equally to the 
superconducting pairing, giving a constant supercon- 
ducting gap A. The difference of A on different sheets 
of the Fermi surface, or multiband superconductivity 
(MBSC), considered theoretically back in the 50s [1], 
has emerged recently as a possible explanation for the 
unusual properties of MgB2 [2] . 

Based on the observation of a sizable difference of A 
on two Fermi surface sheets by angle resolved photoe- 
mission spectroscopy, it has been suggested that the 
layered superconductor NbSe2 is also host to MBSC [3] . 
However, these surface-sensitive measurements were 
performed only down to 5.3 K, close to T c =7.0 K. We 
present evidence of bulk MBSC at low temperatures in 
this compound. 

The thermal conductivity tt of pure samples of NbSe2 
(residual resistivity po=3 /J.Q cm) was measured upon 
warming in a magnetic field (H \\ c axis) [4]. In Fig. 1, 
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k is presented as n/T vs T 2 , enabling a separation of 
the electronic contribution, jco ~ T, and the phononic 
contribution, k 3 ~ T 3 . In zero field, Ko/T goes to zero 
in the T=0 limit as expected for s-wave superconduc- 
tors, thus ruling out the possibility of nodes in the gap 
at any point on the Fermi surface. With the applica- 
tion of a magnetic field, tto/T rapidly increases (Fig. 2, 
main panel) up to H C 2, following closely the increase 
of the electronic specific heat 7 [5]. This behavior is 
very different from that expected for conventional su- 
perconductors. There, electronic quasiparticles remain 
localized inside the vortex cores, such that the increase 
in quasiparticle density (seen as an increase of 7 with 
H) is not followed by an increase of «o. This is shown 
for VaSi in the upper panel of Fig. 2 [6,7]. 

A closer examination of the field dependence of kq 
and of 7 at fields close to H c \ is consistent with the 
presence of localized states, another indication of a 
gap without nodes. Indeed, in a limited field range, kq 
shows an activated increase in field whereas the specific 
heat increases rapidly above H c i (lower inset in Fig. 2). 

In Fig. 3, we show a more detailed comparison of k 
and 7 by plotting the ratio «o/7- In a rough sense, it 
represents the "mobility" of quasiparticle excitations 
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Fig. 1. k/T vs T 2 for NbSo 2 . The line shows the T -> 
extrapolation which yields the electronic contribution kq/T. In 
zero field, kq/T— 0, as expected for a superconductor without 
nodes in the gap anywhere on the Fermi surface. 
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Fig. 3. Field dependence of (main panel). The inset shows 

the field dependence of thermal conductivity for NbSe 2 (circles) 
and VsSi (triangles), plotted as a function of (H — H c \)/H* 
where H* = H c2 /9 for NbSe 2 and H* = H c2 for V 3 Si. 

tex cores observed by muon spin relaxation in NbSe2 
[8]. Moreover, the superconducting coherence length is 
is related to A via £ ~ vf/A, where vf is the Fermi 
velocity. In NbSe2, vf is approximately the same for 
different sheets of the Fermi surface, such that the ratio 
~ 3 gives a ratio of associated superconducting 
gaps A/ A* ~ 3. This value is consistent with previous 
heat capacity and tunneling results [5,9,10]. 

In conclusion, we have identified the anomalous evo- 
lution of thermal conductivity in NbSe2 versus mag- 
netic field with the existence of two length scales in 
the superconducting state. This finds a natural expla- 
nation in a model of multi-band superconductivity, as- 
suming a ratio of larger gap to smaller gap of about 3. 
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Fig. 2. Field dependence of thermal conductivity (circles for 
NbSc2 and empty triangles for X^Si) and specific heat (filled 
triangles for NbSc 2 and squares for VsSi) in NbSe 2 (main panel 
and lower inset) and VaSi (upper inset). The filled circles come 
from a field sweep at 100 mK. Both quantities are normalized 
to the values in the normal state. 

in the vortex state. It is clear that the field dependence 
of no /7 shows a strong change of slope at H* ~ H C 2 /9, 
supporting the existence of a field scale which plays the 
same role for quasiparticle localization as H C 2 does in 
standard s-wave superconductors. In the inset of Fig. 3, 
we show that the low field behavior of NbSe2 can be 
matched to that of VsSi by using H* ~ H C 2/9. 

While £ is associated with the upper critical field 
(£ 2 = <1>u/27rff C 2), a second length scale £* must be 
associated with H* . The crossover between £; and £;* 
with H can naturally explain the shrinking of the vor- 
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